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Abstract 

Rare B decays provide an opportunity to probe for new physics beyond the Stan- 
dard Model. In this paper, we propose to measure the tau polarization in the inclusive 
decay B — > X s t + t~ and discuss how it can be used, in conjunction with other ob- 
servables, to completely determine the parameters of the flavor-changing low-energy 
effective Hamiltonian. Both the Standard Model and several new physics scenarios are 
examined. This process has a large enough branching fraction, ~ few x 1CP 7 , such 
that sufficient statistics will eventually be provided by the B-Factories currently under 
construction. 
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The recent first observation |]J of the inclusive and exclusive radiative decays B — > X s 7 
and B — > K*j have placed the study of rare B decays on a new footing. These flavor 
changing neutral current (FCNC) transitions provide fertile testing ground for the Standard 
Model (SM) and offer a complementary strategy in the search for new physics by probing 
the indirect effects of new particles and interactions in higher order processes. In particular, 
the probing of loop induced couplings can provide a means of testing the detailed structure 
of the SM at the level of radiative corrections where the Glashow-Iliopoulos-Maiani (GIM) 
cancellations are important. These first measurements have restricted the magnitude of 
the electromagnetic penguin transition, resulting in bounds on the value J2| of the ratio 
of Cabibbo-Kobayashi-Maskawa (CKM) weak mixing matrix elements | Vf S / V^i>| , as well as 
providing powerful constraints on new physics H which in some classes of models complement 
or surpass the present bounds obtainable from direct collider searches. 

The study of rare B decays can be continued with the analysis of the higher order 
process B — > X s £ + £~ . The experimental situation for these decays is very promising^], with 
e + e~ and hadron colliders closing in on the observation of exclusive modes with I = e and 
\i final states, respectively. These transitions proceed via electromagnetic and Z penguin 
as well as W box diagrams, and hence can probe different coupling structures than the 
pure electromagnetic process B — > X s 7. Investigation of this decay mode offers exciting 
possibilities as various kinematic distributions associated with the final state lepton pair, such 
as the lepton pair invariant mass spectrum and the lepton pair forward-backward asymmetry, 
can also be measured in addition to the total rate. These distributions are essential in 
separating the short distance FCNC processes from the contributing long range physics^. 

In particular, it has been shown || [7], [8| that the lepton pair forward-backward asymmetry is 
sizable for large values of the top-quark mass and is highly sensitive to contributions from new 
physics. Ali et al. have proposed a program to use these distributions, as well as the total 
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rate for B — > X s 7, to determine the sign and magnitude of each class of short distance FCNC 
contribution in a model independent fashion. Here, we propose a new observable, the tau 
polarization asymmetry for the decay B — > X s t + t~ . We will show that this asymmetry also 
has a large value for top-quarks in the mass range observed at the Tevatron, and will be 
measurable with the high statistics available at the B- Factories presently under construction. 
The tau polarization asymmetry furnishes one more piece of available information for the 
study of rare B decays. Together with the remaining kinematic distributions mentioned 
above, the polarization asymmetry (and the M TT spectrum) will then provide a complete 
arsenal for a stringent test of the SM. 



The transition rate for B —> X s £ + £ , including QCD corrections [10|, is computed via 
an operator product expansion based on the the effective Hamiltonian, 
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H eff = -£v tb V t * s £ dWOiiii) , (1) 

which is evolved from the electroweak scale down to \i ~ m b by the Renormalization 
Group Equations. Here Vij represents the relevant CKM factors, and the 0\ are a com- 
plete set of renormalized dimension five and six operators involving light fields which govern 
b — ► s transitions. This basis (involving left-handed fields only) consists of six four-quark 
operators Ci_6, the electro- and chromo- magnetic operators respectively denoted as 07,8, 
09 ~ es£ Q ,7 A1 &£ Q ,£yY, and O w ~ es Lai ^LcJ-l^l^ '• F° r B — > X s £ + £~ , this effective Hamilto- 
nian leads to the matrix element (neglecting the strange quark mass) 



m = y^v tb v t : 

7T 



(2) 



where q 2 is the momentum transferred to the lepton pair. The Wilson coefficients Cj of the 
b — > s operators are evaluated perturbatively at the electroweak scale where the matching 
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conditions are imposed and are then evolved down to the renormalization scale //. CV_io(Mpp) 
are given by the Inami-Lim functions[|TT||, Cz(Mw) — — 1, and Ci^-q(Mw) = 0. The ex- 
pressions for the QCD-renormalized coefficients Ci(n) are given explicitly in Refs. |7], |10 |. 

The effective coefficient of Og is defined by Cg (/i) = Cg(/i) + ^(a*, <? 2 ) where the function 
Y contains the contributions from the one-loop matrix element of the four-quark operators 
and can be found in Refs. [0, 10 . We note that Y(n,q 2 ) contains both real and imaginary 
contributions (the imaginary piece arises when the c-quarks in the loop are on-shell). The 
differential branching fraction for B — > X s t + t~ is then 
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with all Wilson coefficients evaluated at fi ~ m^, s = q 2 /ml, x = m 2 /m 2 , z = m c /mb, 
and /(^) and represent the phase space and QCD corrections |12[, respectively, to the 
semi-leptonic rate. This agrees with the literature in the zero lepton mass limit. The 
differential branching fraction is scaled to that of the semi-leptonic decay B — > X£u to remove 
the uncertainties associated with the overall factor of ml and to reduce the ambiguities 
involved with the imprecisely determined CKM factors. It is well known that there are large 
uncertainties (up to ±30%) associated with the values of the coefficients CV^/x) due to the 
renormalization scale dependence of the QCD corrections at leading-logarithmic order, as 
well as from the scale parameter in a s . However, this dependence is expected to be reduced 



at the next-to- leading order. This has recently been demonstrated by Buras and Munz[10 



for the case of Cg(/x), which was found to deviate by only ±8% as the renormalization scale 
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\i and the QCD scale parameter Aj^ were varied within their full range of values. The 
situation for C<?(ji) differs, however, as only partial NLO calculations exist. These partial 
calculations do exhibit a reduced /i dependence, and we eagerly await the completion of the 
NLO computations in this case. 

B — ► X s £ + £~ also receives large long distance contributions from the tree-level pro- 
cesses B — > K^ipO followed by ip^ — > £ + £~. These pole contributions are incorporated 
into the lepton pair invariant mass spectrum following the prescription in Ref. f|, where 
both on- and off-shell vector mesons are considered by employing a Breit-Wigner form for 
the resonance propagator. This produces an additional contribution to Cg of the form 

-37T Mypyi - i + n m 

2-^ (S _ A#2 / /rri 2\ i „T_. A/T-. /™2 ' ^ > 



a 2 m 2 



,{s-M 2 Jml)+iT Vi M v Jml 



b Vi=i>,i> 

The relative sign between the short and long distance terms was once a source of controversy, 



but can be explicitly determined via the analyses presented in Ref. |13]. The resulting 
differential branching fraction for B —>■ X s £ + £~ , with and without the long distance resonance 
contributions, is presented in Fig. la for both £ = e and r, taking mt = 180 GeV, mj = 
4.87 GeV, and z = 0.316. We see that the pole contributions clearly overwhelm the branching 
fraction near the ip and ip' peaks, and that there is significant interference between the 
dispersive part of the resonance and the short distance contributions. However, suitable £ + £~ 
invariant mass cuts can eliminate the resonance contributions, and observations away from 
these peaks cleanly separate out the short distance physics. This divides the spectrum into 
two distinct regions^, (i) low-dilepton mass, Ax < s < M^/m 2 — 5, and (ii) high-dilepton 
mass, M 2 i/ml + 5 < s < s max , where 5 is to be matched to an experimental cut. The 
integrated branching fractions (without the pole contributions) for i — e,fi,r are presented 
in Table 1 for both the total and high dilepton mass regions of s. We note that the branching 
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Table 1: Integrated branching fractions for B — > X s £ + £ for the total and high dilepton 
mass regions. 

fraction for B — > X s t + t~ is comparable to that for £ — e, \i in the clean s region above the 
ip' resonance. The exclusive decay B — > Kt + t~ has been computed via heavy meson chiral 
perturbation theory by Du et al. \\14\\ , where the exclusive branching fraction was found to be 
~ 50 — 60% of the inclusive; this places B(B — > Kt + t~) in the range ~ 2 x 10~ 7 . Of course, 
calculations of exclusive decay rates are inherently model dependent ||15||, implying that some 
degree of uncertainty is associated with this result. However, chiral perturbation theory 
is known to be reliable at energy scales smaller than the typical scale of chiral symmetry 
breaking, Acsb — 47r/ 7r /\/2. In B — > Kt + t~, the maximum energy of the K meson in 
the B rest frame is (m B + m 2 K — Am 2 ) /2ms ~ 1.5 GeV, which places most of the available 
phase space at or comfortably below the scale Acsb- We thus expect this method to give a 
reasonable estimate of the exclusive rate. 



The tau polarization asymmetry is defined as 

dB\=-i — dB\=+\ 



P T (s) 



dB x= _! + dB 



A=+l 



(5) 



where dB represents the differential B —>■ X s t + t~ branching fraction. The spin projection 
operator is represented as (1 + 75 $/2, with the normalized dot product being defined as 
s • p = A = ±1 with the — (+) sign corresponding to the case where the spin polarization is 
anti-parallel (parallel) to the direction of the r~ momentum. This corresponds to the usual 
definition of a polarization asymmetry, given in terms of couplings, i.e., (L — R)/(L + R), 
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in the massless case. We note that, of course (and unfortunately), the polarization of final 
state massless leptons cannot be determined in a collider environment. For the process 
B — > X s t + t~ this asymmetry is then calculated to be 



-2 [1 - Ax/s} 1/2 C 10 \RzCl fS {\ + 2s) + 6C 7 
PJs) = — 



D 



(6) 



where D is given by the expression in the curly brackets in Eq. (3). The tau polariza- 
tion asymmetry is displayed as a function of s = q 2 /m 2 in Fig. lb, with and without the 
long distance resonance contributions, and taking m t = 180 GeV. We see that the asym- 
metry vanishes at threshold and grows with increasing s. The value of the total integrated 
asymmetry (i.e., averaged over the high dilepton mass region, s > 0.6) is —0.484. The ex- 
perimentally relevant number of events required to measure an asymmetry a at the na level 
is iV = n 2 /Ba 2 , and is given here by N = n 2 /(4.3 x 10~ 7 )(-0.484) 2 = (n 2 )9.9 x 10 6 for the 
inclusive decay. The exclusive case of B — * Kt + t~ would then yield N ~ (n 2 )2.1 x 10 7 . This 
result demonstrates that P T should be accessible, after several years of running (even when 
r identification efficiencies are taken into account), at the B- Factories under construction. 



The formalism for determining the polarization of r's has been extensively studied fllqj 
for many years (in fact, even before the tau was discovered!). More recently, polarization 
measurements of final state r leptons have been proposed as a useful tool in discerning physics 
beyond the SM in a variety of processes. Some examples include, determining the transverse 
and longitudinal r polarization in B,Af, — > XtvW7§, t polarization asymmetry at hadron 
colliders as a probe of new Z' couplings |TB|, using r polarization to enhance charged Higgs 
boson searches at hadron colliders 0], probing neutralino mixing through r polarization in 



scalar r decay J20|, and in searching for CP violation in the leptonic sector [2]J . These mea- 



surements can take place as information on the tau's polarization state is carried to its decay 
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products. In particular, the momentum distribution of the decay products A (in r — > Au T , 
where A = eu e , \iv^ nn, p, ai, ...) has a large dependence on the spin polarization state of the 
parent r lepton. This dependence is sufficiently striking, such that the tau's helicity may be 
established from a relatively low number of events[IB]. The formalism developed for deter- 



mining the tau's spin polarization on the Z resonance can be applied in this case in the t + t~ 
invariant mass rest frame, except that one cannot take the collinear limit E T ^> m T . The 
resulting decay distributions cannot, however, be measured in this rest frame due to the two 
undetected neutrinos, and one is forced to transform to the laboratory frame to implement 
this procedure. As an example of how accurately the tau's polarization can be determined, 
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we note that the four LEP detectors have made separate polarization measurements |f22 
each of the r decay modes ei/u, pji/u, 7r(K)v, pi/, and a\V (in addition, DELPHI has used an 
inclusive one-prong hadronic analysis). These modes account for ~ 80% of all r decays. The 
tau polarization has then been determined with an overall error of 10 — 15% per experi- 
ment. If B-factory experiments can achieve similar results (and there is no reason to believe 
otherwise) then they will have sufficient statistics to measure the asymmetry P T . 

We now explore the sensitivity of P T to new physics. We first investigate the influ- 
ence of a change in sign of the short distance contributions to CV_io (holding the magni- 
tudes constant). The results are shown in Fig. 2a, where the dashed, dash-dotted, dot- 
ted, solid, and long-dashed curves represent the polarization asymmetry with C W (M W ) — > 
-C 10 (M W ),C 9 , 10 (M W ) -> -C 9jl0 (M w ),C 9 (M w ) -> -Cg(Mw), the SM, and C 7 , 8 (M W ) -> 
—Cj t s(Mw), respectively. We see that there is large sensitivity to any combination of sign 
changes in C 9> iq(Mw), but little variation to a sign change in the electro- and chromo- 
magnetic operator coefficients. This is due to the fact that the operators 0g,io dominate 
the decay in the high s region. We next examine P T in two-Higgs-Doublet models of type 
II, where a charged Higgs boson participates in the decay via virtual exchange in the 7, Z 



penguin and box diagrams. The modifications to the Wilson coefficients in this model are 



given in Deshpande et aZ. [ 23 1 . The resulting tau polarization asymmetry (with s = 0.7 
and m t = 180 GeV) for various values of the charged Higgs mass is presented in Fig. 2b 
as a function of tan/3 = V2/V1, the ratio of vacuum expectation values of the two dou- 
blets. We see that the effect of the H ± is negligible for values of the parameters which 
are consistent with the present constraints from B — > -XyyjJI], |4j> i.e., tan/3 ;> 1 and 
m H ± ;> 240 GeV. Finally, we study the effects of anomalous trilinear gauge boson cou- 
plings in B — > X s £ + £~. The dependence of the Ci(M w ) on these anomalous couplings can 



be found in Ref. |25|]. Figure 2c displays the deviation of P T (for s = 0.7 and m t = 180 GeV) 
with non- vanishing values of the anomalous magnetic dipole and electric quadrupole WW^f 
coupling parameters, Ak 7 and A 7 , respectively, and of the parameter which governs the 
term ig^e tJlUXp {Wj l d\W u — W v d\W^)Z p in the WWZ Lagrangian. For the anomalous cou- 
pling parameters considered here, B —* X s £ + £~ naturally avoids the problem of introducing 
cutoffs to regulate the divergent loop integrals due to cancellations provided by the GIM 
mechanism^]. As expected, we find little sensitivity to modifications in Cj$(My/) from 
anomalous WW'j couplings, but a large variation due to the influence of anomalous WWZ 
vertices in Cq )W (M w ). 

In order to ascertain how much quantitative information is obtainable on the values of 
the Wilson coefficients C^g 10 from the various kinematic distributions, we perform a Monte 
Carlo study. For illustration purposes, we will examine the case where the SM situation 
is realized, i.e., we assume that there is no new physics contributing to these decays. For 
the process B — > X s l + £~ , we consider the M1+1- distribution and the lepton pair forward- 
backward asymmetry for £ = e, /1, and r, as well as the tau polarization asymmetry. 



We take the lepton pair invariant mass spectrum and divide it into 9 bins. These bins are 
distributed as follows: 6 bins of equal size, As = 0.05, are taken in the low dilepton mass 



region below the J ftp resonance, 0.02 < s < 0.32 (where we have also cut out the region 
near zero due to the photon pole), and 3 bins in the high dilepton mass region above the ijj' 
pole, which are taken to be 0.6 < s < 0.7, 0.7 < s < 0.8, and 0.8 < s < 1.0. The number of 
events per bin is given by 



N hi n = C 



/- rf?*' (7) 



and the average value of the asymmetries in each bin is then 

(A) - ^ /^ maX A ^ ' (8) 

-^bin Jsmin ds 

We also include in our study the inclusive decay B — > X s j, which is directly proportional to 
\C 7 (fx)\ 2 . For this case we only consider the total rate. Next we generate "data" (assuming 
the SM is correct) for an integrated luminosity of 5 x 10 8 BB pairs; this corresponds to the 
expected total luminosity after several years of running at future B-factories. The "data" is 
then statistically fluctuated by a normalized Gaussian distributed random number procedure. 

The statistical errors are taken to be 5N = y/N and 5A = \J (1 — A 2 )/N. We include 

statistical errors only for the decay B — > X s £ + £~ , as even for this large value of integrated 
luminosity we expect the errors in each bin to be statistically dominated. The situation differs 
for B — > X s 7, however, as the statistical precision will far exceed the possible systematic (and 
theoretical) accuracy. Hence, in this case we assume a flat 10% error in the measurement 
of the branching fraction. We then perform a three dimensional x 2 fit to the coefficients 
CV,9,io(a0 from the "data" according to the usual prescription 

/nobs _ nSM\ 2 

where Q° bs,SM , 5Qi represent the "data", the SM expectations, and the error for each observ- 
able quantity Qi. The resulting 95% C.L. allowed regions as projected onto the Cg(/i)— Cio(/x) 



and Cr(fJ-) — Cio(n) planes are presented in Fig. 3a and b, respectively. In these figures, the 
point 'S' labels the SM expectations (assuming m t = 180 GeV), and the diamond represents 
the best fit value which has a total x 2 — 24.6/25dof. We see that the determined ranges for 
the coefficients encompasses their SM values. The 95% C.L. ranges for the coefficients are 
found to be C 7 (n) = 0.32081^11 , C g (fj) = -2.3001^1, and C 10 (fj) = 4.834±g;f^, corre- 
sponding to a 7.5% , 20% and 10% determination of C 7 , C 9 , and C w , respectively. Clearly, 
the values extracted for C$ and Cio are highly correlated, but this is not the case for C 7 and 
Ci . If we take C 7 {ji) to have opposite sign, i.e., C 7 (/i) = — \C 7 (fi)\, we find that the fit is 
quite poor with the best fit yielding x 2 — 539.1/25dof. 

In conclusion, we have shown that measurement of the r polarization in B — > X s t + t~ 
is highly sensitive to new physics and hence provides a powerful probe of the SM. Together, 
measurement of the polarization asymmetry and the remaining kinematic distributions as- 
sociated with B — > X s £ + £~, will provide enough information to completely determine the 
parameters of the FCNC effective Hamiltonian. We find that the values of the polarization 
can be precisely determined with the large data samples that will be available at the B- 
Factories presently under construction. We eagerly await the completion of these machines! 

ACKNOWLEDGEMENTS 

The author thanks T.G. Rizzo for invaluable discussions, and the Phenomenology 
Institute at the University of Wisconsin for their hospitality while this work was completed. 



10 



References 



[1] M.S. Alam et al, (CLEO Collaboration), Phys. Rev. Lett. 74, 2885 (1995); R. Ammar 
et al, (CLEO Collaboration), Phys. Rev. Lett. 71, 674 (1993). 

[2] A. Ali and C. Greub, Z. Phys. C60, 433 (1993). 

[3] For a review see, J.L. Hewett, in Proceedings of the 21st Annual SLAC Summer Institute, 
ed. L. DePorcel and C. Dunwoodie (SLAC-Report-444, Stanford, CA, 1994). 

[4] R. Balest et al, (CLEO Collaboration), to appear in Proceedings of the 1994 Inter- 
national Conference on High Energy Physics, Glasgow, Scotland, CLEO-CONF-94-4 
(1994); C. Anway-Wiese, (CDF Collaboration), in Proceedings of the 8th Meeting of the 
Division of Particles and Fields of the American Physical Society, Albuquerque, New 
Mexico, 1994, ed. S. Seidel (World Scientific, Singapore, 1995). 

[5] N.G. Deshpande, J. Trampetic, and K. Panrose, Phys. Rev. D39, 1461 (1989); C.S. 
Lim, T. Morozumi, and A.I. Sanda, Phys. Lett. B218, 343 (1989). 

[6] A. Ali, T. Mannel, and T. Morozumi, Phys. Lett. B273, 505 (1991). 

[7] A. Ali, G.F. Giudice, and T. Mannel, Z. Phys. C67, 417 (1995). 

[8] C. Greub, A. Ioannissian, and D. Wyler, Phys. Lett. B346, 149 (1995); D. Liu, Univ. 
of Tasmania Report (1995); G. Burdman, Fermilab Report Pub-95/113-T (1995). 

[9] F. Abe et al, (CDF Collaboration), Phys. Rev. Lett. 74, 2626 (1995); S. Abachi et al., 
(DO Collaboration), ibid., 74, 2632 (1995). 



11 



[10] N.G. Deshpande and J. Trampetic, Phys. Rev. D60, 2583 (1988); B. Grinstein, M.J. 
Savage, and M.B. Wise, Nucl. Phys. B319, 271 (1989); A.J. Buras and M. Miinz, 
Phys. Rev. D52, 186 (1995). 

[11] T. Inami and C.S. Lim, Prog. Theor. Phys. 65, 297 (1981). 

[12] N. Cabibbo and L. Maiani, Phys. Lett. 79B, 109 (1978). 

[13] P.J. O'Donnell and H.K. Tung, Phys. Rev. D43, R2067 (1991); A.I. Vainshtein, V.I. 
Zakharov, L.B. Okun, and M.A. Shifman, Yad. Fiz. 24, 820 (1976) [Sov. J. Nucl. Phys. 
24, 427 (1976)]; N.G. Deshpande, X.-G. He, and J. Trampetic, Univ. of Oregon Report 
OITS-564 (1995). 

[14] D. Du, C. Liu, and D. Zhang, Phys. Lett. B317, 179 (1993). 

[15] For a review, see R. Grigjanis, P.J. O'Donnell, M. Sutherland, and H. Navelet, 
Phys. Rep. 228, 93 (1993). 

[16] Y.S. Tsai, Phys. Rev. D4, 2821 (1971); T. Hagiwara, S.-Y. Pi, and A. Sanda, Ann. 
Phys. (NY.) 106, 134 (1977); J. Kiihn and F. Wagner Nucl. Phys. B236, 16 (1984); P. 
Aurenche and R. Kinnunen, Z. Phys. C28, 261 (1985); K. Hagiwara, A.D. Martin, and 
D. Zeppenfeld, Phys. Lett. B235, 198 (1990); A. Rouge, Z. Phys. C48, 75 (1990). 

[17] A.F. Falk, Z. Ligeti, M. Neubert, Y. Nir, Phys. Lett. B326, 145 (1994); Y. Grossman 
and Z. Ligeti, Phys. Lett. B347, 399 (1995); R. Garisto, Phys. Rev. D51, 1107 (1995); 
M. Gremm, G. K6pp, and L.M. Sehgal, Phys. Rev. D52, 1588 (1995). 

[18] J.D. Anderson, M.H. Austern, and R.N. Cahn, Phys. Rev. Lett. 69, 25 (1992), and 
Phys. Rev. D46, 290 (1992). 



12 



[19] B.K. Bullock, K. Hagiwara, and A.D. Martin, Phys. Rev. Lett. 67, 3055 (1989); 
Nucl. Phys. B395, 499 (1993); D.P Roy, Phys. Lett. B277, 183 (1992); S. Raychaudhuri 
and D.P. Roy, Phys. Rev. D52, 1556 (1995). 

[20] M.N. Nojiri, Phys. Rev. D51, 6281 (1995). 

[21] Y.S. Tsai, Phys. Rev. D51, 3172 (1995). 

[22] D. Buskulic et ai, (ALEPH Collaboration), CERN Report CERN-PPE/95-023 (1995), 
and Z. Phys. C59, 369 (1993); P. Abreu et ai, (DELPHI Collaboration), Z. Phys. C67, 
183 (1995); M. Acciarri et al, (L3 Collaboration), Phys. Lett. B341, 245 (1994); R. 
Akers et al, (OPAL Collaboration), Z. Phys. C65, 1 (1995). 

[23] N.G. Deshpande, K. Panrose, and J. Trampetic, Phys. Lett. B308, 322 (1993). 

[24] J.L. Hewett, Phys. Rev. Lett. 70, 1045 (1993); V. Barger, M. Berger, and R.J.N. Phillips, 
Phys. Rev. Lett. 70, 1368 (1993). 

[25] S.-P. Chia, Phys. Lett. B240, 465 (1990); K.A. Peterson, Phys. Lett. B282, 207 (1992); 
T.G. Rizzo, Phys. Lett. B315, 471 (1993); G. Baillie, Z. Phys. C61, 667 (1994). 

[26] We note that the expressions given for the lepton pair forward-backward asymmetry in 
Refs. 6 and 7 do not agree; we have found agreement with that of Ref. 6. In the case of 
massive leptons, i.e., i = r, one must include an additional overall factor of [1 — Ax/s] 1 ^ 2 
in the asymmetry. 



13 



X) 



I 

+ 

w 

X 

t 

PQ 

> — ' 

m 




s = q /m b 



0.4 



0.2 



0.0 



-0.2 



-0.4 



-0.6 



-0.8 



n — i — i — i — I — i — r 



(b) 



_i i i i i i i_ 



_i i i i i i i i_ 



0.5 



0.6 



0.7 



0.8 



0.9 



2/2 

s = q /m 



Figure 1: (a) Differential branching fraction and (b) tau polarization asymmetry as a function of 
s for £ = t (solid and dashed curves) and I = e (dotted and dash-dotted curves) , with and without 
the long distance resonance contribution. 
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Figure 2: Tau polarization asymmetry (a) with changes in the sign of the Wilson coefficients at 
the electroweak scale, corresponding to Cio, Cg^o, Cg, SM Cj t s from top to bottom; (b) in two- 
Higgs-doublet models as a function of tan/3 with m H ± = 50,100,250,500 corresponding to the 
solid, dashed, dotted, and dash-dotted curves, respectively. The SM value is denoted by the solid 
horizontal line, (c) with anomalous couplings WWj and WW Z couplings as described in the text. 
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Figure 3: 95% C.L. contour in the (a) Cg — C±o, (b) C-j — C10 plane from the fit procedure described 
in the text. 'S' labels the SM prediction and the diamond represents the best fit values. 
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